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Abstract 

Nonlinear conservation laws driven by Levy processes have solutions which, in 
the case of supercritical nonlinearities, have an asymptotic behavior dictated by the 
solutions of the linearized equations. Thus the explicit representation of the latter is 
of interest in the nonlinear theory. In this paper we concentrate on the case where the 
driving Levy process is a multiscale stable (anomalous) diffusion, which corresponds 
to the case of multifractal conservation laws considered in [3 H [31E]. The explicit 
representations, building on the previous work on single-scale problems (see, e.g.,[5]), 
are developed in terms of the special functions (such as Meijer G functions), and are 
amenable to direct numerical evaluations of relevant probabilities. 


1 Introduction 

Mathematical conservation laws are integro-differential evolution equations, such as Navier- 
Stokes and Burgers equations, expressing the physical principles of conservation of mass, 
energy, momentum, enstrophy, etc., in different dynamical situations. With this paper 
we initiate a program of investigation of explicit representations for asymptotic behavior 
of solutions of conservation laws driven by multiscale, (ai,..., afc)-stable. Levy processes 
(multifractal anomalous diffusions). Since the asymptotic behavior of such conservation 
laws is determined, in some cases, by their linearized versions, the starting point here is to 
obtain exact representation, via known special functions such as Meijer G functions, of the 
solutions of linear multiscale evolution equations, that is, for the PDFs of the multiscale 
Levy processes themselves. 

The idea is to produce a framework that permits a straightforward calculation of 
probabilities related to the mutiscale diffusions using a symbolic manipulation platform 
such as Mathematica, and a fairly standard set of special functions that have been in use 
in this area for a long time. 

The plan of the paper is as follows. We begin, in Section 2, with a review of the known 
general results from mniias! on the asymptotics of solutions of multifractal conservation 
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laws, and apply them to the case of general asymmetric two sided multiscale diffusion. 
In the case of supercritical nonlinearity asymptotics is dictated by the linear part of the 
equation so, in Section 3, we produce an exact representation of solutions of linearized 
equations in the general asymmetric case; simpler representation are then deduced in the 
symmetric case. In Section 4, anticipating our future needs to obtain explicit solutions 
for equations describing subdiffusive anomalous diffusions, where the time is also ’’fractal” 
[6], we obtain an explicit representation for totally asymmetric a-stable diffusions, for 
0 < a < 1. Conclusions, as well as a discussion of the relevant moment problem, can be 
found in Section 5. 

In the remainder of this section we establish the notation and provide basic definitions 
of the integral transforms and special functions we are going to work with. We start with 
the Fourier transform of an integrable function /(x) defined for x G M, and real oj [7], 

/ CX> 

e^‘^^f{x)dx, ( 1 . 1 ) 

-OO 

1 

fix) = -F-i[/(a;);x] = ^j_ e--V(^)da;. (1.2) 


For a function /(x) = 0, for x < 0, such that e ^^f{x) is integrable on the positive half-line 
for some fixed number c > 0, the Laplace transform 


roo 

f*ip) = ^[fix);p] = / e-P^fix)dx, 

Jo 

1 rc^ioo 

fix) = c-^[npy,x] = — / epy*ip)dp, 

J c— ico 


(1.3) 

(1.4) 


where p = c + iuj. There is an obvious relationship between the Fourier transform of 
/i(x) = e““/ 2 (x) and the Laplace transform of / 2 (x), see e.g., [7] for more information. 
Finally, the Mellin transform of /(x) is here defined as follows: 

poo 

fis) = M[fix)\s\= j x^~^fix)dx, (1.5) 

xo 

fix) = M~^[fis)-x] = ^ (1-6) 

where s is a complex variable [7]. The contour of integration L is determined by the 
domain of analyticity of f 2 is) and, usually, it is an infinite strip parallel to the imaginary 
axis. The conditions under which the integrals in Eqs. (jl.ip - (jl.6p converge can be found 

in [7]. 

The Meijer G function will play a pivotal role in what follows. It is defined as the 
inverse Mellin transform of products and ratios of the classical Euler’s gamma functions. 
More precisely, see [Hi, 


/^m,n 

^p,q 




'p \ _ 


= M- 


in%m+i r(i - Bj - s) U%n+1 +«) ’ 


; 2 


(1.7) 


where empty products in Eq. (jl.Tjl are taken to be equal to 1. Eq. (II.7h holds under the 
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following assumptions: 


z/0, 0 <m < q, 0 < n < p, 

AjeC, j = l,...,p; BjeC, j = l,...,q. (1.8) 


A description of the integration contours in Eq. (HZD, and the general properties and 
special cases of the Meijer G functions can be found in [5]. If the integral in Eq. Q 
converges and if no confluent poles appear among r(l — Aj — s) or r(l — Bj — s), then 
the Meijer G function can be expressed as a finite sum of the generalized hypergeometric 
function, see formulas (8.2.2.3) and (8.2.2.4) on p. 520 of [8]. Recall that a generalized hy¬ 
pergeometric function can be represented in terms of the following series, see Eq. (7.2.3.1) 
on p. 368 of [8]: 


/ ai, . . . , Qp \ x” Wj=l{0‘j)n 


(1.9) 


where the upper and lower lists of parameters are denoted by (op) and (6^), respectively, 
and {a)n = r(a -|- n)/r(a) is the Pochhammer symbol. 

To conclude the introduction we find it convenient to introduce the special notation 
for a specihc uniform partition of the unit interval. 


A(n,a) = {^,^,...,^}. (1.10) 

For later reference, we also quote the Euler’s reflection formula, 

r{z)T{i-z) = -^, (1.11) 

sm(7rz) 

see Eq. (8.334.3) on p. 896 in |10] . and the Gauss-Legendre multiplication formula. 


r(na) 


n—1 

( 27 r) r 


j=0 



( 1 . 12 ) 


see Eq. (8.335) on p. 896 of [TO] . 


2 Multiscale Levy processes 


In this Section we are turning to a review of infinitesimal generators A of semigroups 
associated with 1-D multiscale (ai,..., Q:fc)-stable Levy processes driving the evolution 
equations of the form. 


du 9 

-+Au+-g{u)^0, 


u{0,x) = no(x). 


( 2 . 1 ) 


where u = u{t,x),t > 0,x G M, uq is an initial condition which will be specified later, and 
5( : M I— )• M is a (nonlinear) function. Such equations are often called fractal, or anomalous 
conservation laws Their asymptotic behavior will be discussed in Section 3. 

The operators A are easiest to describe in terms of their actions in the Fourier do¬ 
main; they are so-called Fourier multiplier operators. Let us begin by recalling the basic 
terminology and establishing the notation. 
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Like any Markov processed, the Levy process, Xt,t > 0, has associated with it a 
semigroup Pt of convolution operator^ acting on a bounded function 0(x) via the formula, 


tcpix) =E^{(j){X{t)) = [ (l){x + y) P{X{t) e dy). 
Jr 


The infinitesimal generator A of such a semigroup is defined by the formula, 

Ph-Po 


— A = lim 


/n>o h 


( 2 . 2 ) 


(2.3) 


and the family of functions,u(t, x) = Ptg{x), interpreted here as probability density 
functions (PDFs), clearly satisfies the (generalized) Fokker-Planck evolution equation. 


dv . 

at = 

because lim{Pt+h - Pt)/h = lim[(P,j - Po)/h]Pt = -APt- 

h^O h^O 

In the case of a general Levy processes Xt, we have the identity. 


(2.4) 




(2.5) 


where P stands for the Fourier transform, and 

V’(w) = logE[e*‘^^i] (2.6) 

is the characteristic exponent of Xi, which is necessarily (see, e.g., [12] )j of the form 

= iy,uj - + [ (e*^"^ - 1 - iujxl\^\^i)A{dx), (2.7) 

where y, £ M, a € M+, and A is a nonegative measure on M, satisfying the conditions 
A({0}) = 0, and Jjjj(l A |a;p)A((ix) < oo. The triplet {y,a,A) is called the characteristic 
triplet of X, y £ M - the drift coefficient, a > 0 - the Gaussian, or diffusion coefficient, 
and A - the Levy measure of Xi. The Levy measure describes the “intensity” of jumps 
of a certain height of a Levy process in a time interval of length 1. Observe that 

P{Ptfi){oj) = e-^^^M{Xt + x) dx^ dy 

_^yLUjXt f e“*‘^^(/)(y) dy = exp(ti/'(u;))J'(/)(a;), 

Jr 

which, in view of Eq. (j2.3|) . indeed implies Eq. (12.5h . 

In the case of the usual Brownian motion the infinitesimal operator A is just the 1-D 
classical Laplacian A (the second derivative operator). Eor the self-similar (single-scale) 
symmetric a-stable process Xt, the infinitesimal generator is the 1-D fractional Laplacian 
_(_A)«/2, 0 < q; < 2, corresponding to the characteristic exponent (Eourier multiplier) 
fi’iuj) = —|a;|". 

^See, e.g., m, for basic information in this area. 

^ That is, Pt+s = PtPs, L s > 0. 
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In what follows we focus our attention on the multiscale (and not necessarily symmet¬ 
ric) Levy processes with the characteristic functions of the form, 

n 

= ( 2 . 8 ) 


where, for each j = 1,2,... , n. 






(2.9) 


The symbol sgn(a;) denotes the sign of the parameter oj. The multiparameter = 

(tti, /3i,.../3n; 71 ,...) 7n) Las to satisfy the following conditions: If 0 < < 1 

then |/3j| < a^, and if 1 < Oj < 2 then \f5j\ < 2 — a^; for all j, we assume that 7 j > 0 . 
Thus the Fourier multiplier describing the infinitesimal generator of Xt is of the form, 

n 

= E -7i|a;r^'ef (2.10) 

i=i 

The generator itself will be denoted For the sake of convenience, and without 

loss of generality, in the remainder of the paper we will assume that 


ai < a2 <■■■< an 


The densities Vaj,i 3 jix,t) appearing in (j2.9l) are unimodal [l3l[Il]. The skewness pa¬ 
rameter, /3j, measures the degree of asymmetry of Vaj,i 3 j{x,t): for /3j = 0 they are just 
the previously mentioned symmetric aj-stable densities with fractional Laplacians as the 
corresponding infinitesimal generators. Moreover, all of those densities are self-similar, 
since, for any a; G M, and t > 0, 


Vaj,0j,'yj {ti x) 


1 



( 2 . 11 ) 


and 

Vaj,/3j,'yji~X,t) = Vci,^^-/3^^'Yj{t,x). (2.12) 

Eq. (12.1211 is a consequence of the identity —co) = satisfied by 

given in Eq. ([2J]). 


3 Asymptotics of solutions of multifractal conservation laws 
with supercritical nonlinearity 

Now, we are ready to state the results about existence, uniqueness, and the asymptotic 
behavior of the solutions of the Cauchy problem for the multifractal conservation laws 
(j2.ip driven by multiscale Levy processes (anomalous diffusions) introduced in Section 
2. The main point here is the observation that the solutions of Eq. m, under certain 
conditions on the generator A and the nonlinearity g, have the large time behavior similar 
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to solutions of Eq. (j2.4l! l^. These results provide the motivation for the work presented in 
the following sections. The physical justification for considering conservation laws driven 
by Levy processes are more numerous than can be cited here, but see, e.g., m for a review 
of the subject. 

The solutions of Eq. (|2.ip have to be understood in some weak sense which opens 
several possibilities presented for example in nisi Em [19]. Motivated by the classical 
Duhamel formula we choose to interpret them as the so-called mild solutions satisfying 
the identity, 

u{t,x) = [e~^-^uo]{x) - j -^e~^^~'^'^-^[g{u)]{T,x)dT. (3.1) 

The basic results are summarized in the following Theorem, where the regularity of 
the solutions of Eq. (13.ip is expressed in terms of the Sobolev space W^’^. 


Theorem 1 fsee [3]J fij Assume that g G and A is the infinitesimal generator 

of a Levy process with the symbol satisfying the condition 


lim sup 

|lj|—)- oo 


V'(w) - 


a; 


< oo 


for some 0 < a < 2, and ifo > 0. 


(3.2) 


Given uq G L^(M) n L°°(M), there exists a unique solution u G C([0, oo); L^(M) n L°°(M)) 
of the problem 

On 

— + Au+ V ■ g{u) = 0, u{x,0) = uo{x). (3.3) 

This solution is regular, u G (^((O, oo); IT^’^(M)) n C^((0, oo);L^(]R)), satisfies the conser¬ 
vation of integral property, f u(x,t)dx = f uo(x)dx, and the contraction property in the 
Lp(M) space, 

||'w(t)||p < ll^ollp) (3-4) 

for each p G [l,oo], and all t > 0. Moreover, the maximum and minimum principles hold, 
that is, 

ess inf tto < u{x,t) < ess suprto, a.e. x,t, (3.5) 

and the comparison principle is valid, which means that if uq < uq G L^(M), then 


u{x,t) < v{x,t) a.e. x, t, and ||n(t) — u(f)||i < ||tto — uo||i. 
(ii) Under the following additional conditions on the symbol of A, 

V'(w) 


, , ibioo) , L’(w) 

0 < hmmf —r— < hmsup -j—j— < oo. 


0 < inf ■ 


LO \UJ 


uj->o |a;|“ |w, 

for some 0 < a < 2, the more precise bound, 

||«(t)||p < 

holds for all 1 < p < oo. Moreover, if uq G L^(M) n then 


2 ’ 


(3.6) 


(3.7) 


(3.8) 


®This is in contrast to the phenomena observed for data of Riemann type (nonintegrable, and nons¬ 
mooth), when shocks are created, see .e.g., |15l 1161 fT7| . 
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with a constant C which depends only on ||mo||i ||tto||p- 

(in) Assume that u is a solution of the Cauchy problem (13.31) with uq G L^(M) n 
and that the symbol if of the generator A satisfies Eqs. (13:^ and (13771) with some 0 < 
a < 2. Furthermore, suppose that the nonlinearity g is supercritical, that is, g £ , and 

limsup^^o lfl'('S)l/|sr < oo? some r > max(a, 1). Then the relation 

lim - e-^-^uoL = 0 (3.9) 

t^oo 

holds for every 1 < p < oo. As usual, denotes the action of the Levy semigroup on 

the function uq, i.e. is a solution of the linear Eq. dMD with the initial data uq. 

On the other hand, the asymptotics of the solution of the linear Cauchy problem 
Eq. (|2.4p is well known: there exists a nonnegative function g G L°°(0, oo) satisfying 
lim g{t) = 0 such that 

£—)-oo 


< f(i"i/p)/"?7(t), (3.10) 

p 

where p^(t) is the kernel of the operator A in Eq. (12.3p . Higher order asymptotics is also 
available [2]. 

The above general results have direct consequences for multifractal conservation laws 
driven by multiscale anomalous diffusions introduced in Section Note the parabolic 
regularization included in the operator A because of the conditions (13.21) . and (|3.7I) . 


_ f A 

e * uo 


- / uo{x) dx ■ pA{t) 


Corollary 2 All the statements of Theorem are valid for the conservation laws 


du d 


(3.11) 


with 


a = Ol < 02 < • • • < On = 2. 

In particular, if u is a solution of the Cauchy problem (j3.1ip withuQ G L^(M)nL°°(M), and 
the nonlinearity g G is supercritical, i.e., limsup^^Q | 5 ((s)|/|s|’' < oo, forr > max(Q;, 1), 
then the relation 


lim 

t—^OO 


— e 


—tA 


'^ollp — 0 


(3.12) 


holds for every 1 < p < oo. Moreover, 





dx ■ PA - (t) 


p 


where p_A - (t) is the kernel of the operator A^ - in Eq. (|3.11|). 


^The particle approximations and the propagation of chaos results for such systems have been studied 
in [20]. 
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To prove Corollary [2] it suffices to show that conditions Eq. (13.21) and Eq. (13.7p are 
satisfied. Indeed, for the mutiscale Levy process with symbol (I2.10p . we have 

M^oo ' 

with aj* = a* where a* = max(ai,..., Also, 

0 < hm —;—;-= a,„ < oo, 

with aj, = a* where a* = min(ai,..., a^) = ai; and 

inf - > 7n > 0. 

cur 


The above results depended on the subcritical behavior 

limsup < oo, for r > max(a, 1), 

s^o |sr 

of the nonlinearity in the conservation laws discussed above. Note that in the classical 
case of the Burgers equation the situation is dramatically different. 

Remark 1. Asymptotics of solutions of the Burgers equation. The first order asymp¬ 
totics of solutions of the Cauchy problem for the Burgers equation 

is described by the relation 

t^^-^/^'^/^\\u{t) — UM{t)\\p ^ 0, as t —)• oo, 

where 

g-3:V(4t) / I r^Vi 

UM{x,t) =iK{M) + - e^l^doj 

is the so-called source solution with the initial condition u(x, 0) = MSq. It is easy to verify 
that this solution is self-similar, i.e., UM(x,t) = t~^^^U(xt~^/^,l'). Thus, the long time 
behavior of solutions of Eq. (|3.13p is genuinely nonlinear, i.e., it is not determined by the 
asymptotics of the linear heat equation. This strongly nonlinear behavior is due to the 
precisely matched balancing influence of the regularizing Laplacian diffusion operator and 
the gradient-steepening quadratic inertial nonlinearity, see [2III221II8]. 

Although not needed explicitly in the remainder of the paper, for the sake of com¬ 
pleteness we are providing below a general result showing how such a matching critical 
nonlinearity exponent for the nonlocal multifractal conservation law yields the solutions 
of (13.3p which behave asymptotically like the self-similar source solutions U of ()3.3p with 
singular initial data MSg- 












Theorem 3 (see m) Let 1 < a < 2, and u be a solution of the Cauchy problem i3.3\) with 
the operator A = (—+ JC, with the perturbation K, being another Levy infinitesimal 
generator whose symbol k fulfills the condition, 


lim 



= 0 , 


(3.14) 


and uq G uo{x) dx = M > 0. Assume that g satisfies the condition 


Then, for each 1 < p < oo. 


r 9 {s) 

lim —— 7 —^ 


G 


(3.15) 


lim t^^-^M/^u{t) - U{t)\L = 0, 

t^CXD 


(3.16) 


where U = Um is the unique solution of the problem 113. 51) with r = a and the initial data 
M 6 q. Moreover, U is of self-similar form U{x,t) = 1), J^dLf{x,l)dx = 

M, and U >Q. 


Thus, analogous to Corollary [2l we also have the following result in the case of multifractal 
conservation laws with critical nonlinearities. Note that, in contrast to Corollary [2l the 
parabolic regularization is not necessary here. 

Corollary 4 All the statements of Theorem\^are valid for the multifractal conservation 
laws 

Ou 0 

'^9{u) = 0 , u{x,0) = Uo{x), (3.17) 

with a = = min(Q;i,..., a^) = ai. 

The verihcation of the condition (j3.14p is immediate. With the symbol of the pertur¬ 
bation /C, 

n 

k{w) = 

i =2 

we do have lim k{uj)/\ujC = 0. Recall that, in view of the convention adopted at the 
beginning of the paper, a* = oi < 02 < • • • On- 

The issue of explicit representations of source solutions of fractal conservation laws 
with critical nonlinearities is obviously more difficult than the problems we are addressing 
in the subsequent sections, but we plan to investigate it in the future. 


4 Explicit representation of the kernels of the two-scale, 
two-sided Levy generators, 0 < ct < 2 

In this section our goal is to find explicit representations for kernels of the infinites¬ 

imal generators A^^. 0 ^^ which dictate the long-time behavior of the nonlinear conservation 
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laws discussed in Section 3. For the sake of simplicity, we present the case when the scal¬ 
ing parameter 7 = the notation is then streamlined to Simply 

stated, we need to find an explicit expression for the Fourier convolution of 
j = 1 , 2 , X € M, and t > 0: 

roc roo 

H{t,x)= - y)dy = Vai,p:,(t,x - y)va2,f32{^,y)dy, (4.1) 


where H{t,x) = H{d; fi-,t,x), a = (ai,a 2 ), and f5 = (/ 3 i,/ 32 ). The basic properties of 
with necessary conditions on aj and fij, are given in Section 2. 

The functions Vaj,Pjit,x), j = 1,2, represent the unimodal probability density func¬ 
tions of two-sided Levy stable distributions Haul], which correspond to one-sided Levy 
stable distributions for 0 < aj < 1 and fdj = —cxj. This case will be discussed in section 
5. The series representation of two-sided Levy stable distributions for 0 < aj < 1, and 
|/3j| < aj, can be found in, e.g., Eq. (5.8.8a) on p. 142 [13], Eq. ( 6 . 8 ) on p. 583 of [Hj . 
and Eq. (4) in [23], whereas, for I < aj < 2, and \/3j\ < 2 — aj, they are described in, 
e.g., Eq. (5.8.8b) on p. 142 of [l3], and Eq. (6.9) on p. 583 of [Hj. Those two different 
types of series expansions were calculated for rational values of parameter aj and /3j, see 
Eqs. (4) and (5) in [5|. We quote some solution which will be used later in the paper: the 
Gaussian distribution 

exp(-lJ) 


V2fl{t,x) = 


2 y/^ 


for a = 2 and /3 = 0, the Levy-Smirnov distribution 


^exp(-ll) 


VI i{t,x) = 

2 ’ 2 2^'kx'^I^ 

= 0 , 


X > 0 

X < 0. 


(4.2) 

(4.3) 

(4.4) 


for a = 1/2 and /3 = — 1 / 2 , and 

(2/t)2/3F(|) ^ /5/6. 4x3 \ , (2/t)4/3r(|) 


30F r(|)^ ^U/3’ 27t2y)+ 90F r(|) ^ ^U/3’ 27t2y) 

/ CO 

g-icjXg-t|(^|3/2exp(-2fsgn(tj)^^ ( 

-00 




for a = 3/2 and /3 = —1/2 |5]. The symbol lEi stands for the hypergeometric function 
introduced in Section 1. 

Let us now hnd the explicit form of H{t,x) given in Eq. (14.11) . Applying the property 
(I2.12P to Eq. (14.ip . we can rewrite H in the form. 


H{t,x) = F ^[h„i,, 9 i(La;)h„ 2 ,/ 32 (Lw);a:] = -x)0(-x) -h iL+(t, x)0(x), (4.6) 


where 


H+{t,x) = H+{ai,/3i,a2,P2;t,x) = 


^ r 

TT Jo 


e-ixu^ exp 


(4.7) 
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and 


= H+{ai,-l3i,a2,-/32;t,x). (4.8) 

The function Q{x) is here the usual Heaviside step function. The matching, at x = 0, 
of these two components is assured by the continuity at the origin of H{t,x), and of all 
of its higher derivatives. Indeed, the continuity of Eq. (j4.6jl at x = 0 can be shown by 
employing Eq. (14.71) as follows: for n = 0,1, 2,... 


^ Jo 

_ }_ f g— cos n)] —ta;“2 cos [^(/32—n)] 

TT Jo 

X COS sin [f (/3i — n)] + sin [f (/32 — ?^)]) du 

= Jh “P (-fexp (-f/ 32 )^^ 

^ Jo 

= '9x-H'-(T-a^)L=o- 


(4.9) 


We would also like to point out that Eqs. (14.6p and (14.71) . in the case ai = a 2 -, and fJi = (32, 
imply the identity, 

H{x,t) = Vai,i3i{x,2t). (4.10) 

In what follows, without loss of generality, we will consider only the case of H+(x,t). 
The function H_(x,t) will be used only when necessary. We assume that, for certain 
values of complex s, the Mellin transform of 77+(x, t) exists and, according to the notation 
introduced in Eq. (jl.5l) . it is denoted by H+{s,t). Thereafter, we substitute Eq. (14.7h into 
77+(s,t), change the order of integration and use the one of Eqs. (2.3.2.13) of |24] . Those 
steps imply that for rational aj and f3j, j = 1,2, such that ai = ^, /3i = a 2 = 

and (32 = where I, k, p, q, a, and b, are integers, we have 


«+(<'>*) = X 77 E 


Mi -1 




j! 


j=0 
^ l+mi 


^ Ml —1 oo 


M M 

j. —+ ( —-Dt 
t M >■> 


l,A(mi,7^ + fj) 

A(Mi,l+j) 
i^_iy+rM, r(,)r[i^ + f(j+rMi)] 


/_ Ml / mi ) \ 

’ V h'h ' j I 


Mtt ^(j + rMi)! ^^+(t-i)(i+rMi) 

X sin [vrai^ - 7r(j +rMi)(u - u^)] 

1 ~ (- 1 )^ t-T-^(M-i)r(s)r(7-* 

^ho~^ r[l-n^ + (u-nf)r]r[u^-(n-uf)r] 


_+ "iri 

M ^ M T 


(4.11) 


(4.12) 

(4.13) 


where m, M, mi, and Mi are as follows: 


m = min(^,-), M = max(^,-), mi = min(/cp,/g), and Mi = msix{kp,lq). 

(4.14) 
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The parameters u, and v, are determined by the equalities, 

u = y, V = -, for ai > a2, and u = -, v = y for ai < a2- (4-15) 

k q q k 

In Eq. ()4.11l) we utilized a series representation of the generalized hypergeometric function 
given in Eq. (II.9p . and the Gauss-Legendre multiplication formula defined in Eq. (I1.12p . 
In Eq. (I4.12P we applied Eq. (11.111) . and also changed the summation index as follows: 
j + rMi —)• r. 

The next step requires inverting the Mellin transform in Eq. (14.131) . To accomplish 
this task we will introduce the new variable of integration, s = (1 — s)/{lp), Eqs. (jl.lip . 
and (|1.12l) . Putting the all of these terms together, we get, for x > 0, 


H+{x,t) 


mi 


y/M 


x{27r 




E 

r=0 


i-ty 


r\ 




^lp-\-u7ni ,mi -\-umi 


/(^Ipypt^i A{lp,0),A{umi,-{v-u^r)) \ 
[xiPmT^ A(mi,fr),A(nmi,-(u-nfr))J 


(4.16) 


The Meijer G functions in Eq. (I4.16P can be expressed, via formulas (8.2.2.3), and (8.2.2.4), 
on p. 520 of [8], in terms of a generalized hypergeometric function. With respect to the 
values of ai and 02 ; we can consider two different cases: 


(A) After applying Eq. (8.2.2.3) on p. 520 of [8] to Eq. (I4.16p . for 0 < ctj < 1, i = 1,2, 
we have 


rr / 1 r(l + Mj + mr) 

H+{x,t) = --}_^ - -sm(r7rn + j7ra) 

^=0 j=o 


^ l+lpP'mi 


j,l+Mj+rm 

1, A{lp, 1 + Mj + mr) 


;(-i) 


miu—mi 






mT^x^P J ’ 


(4.17) 


A(mi,l+j) ,, 1 .^ 

where u, and v, are given in Eq. (j4.15(^1 . Moreover, using the series expansion of 
the function lyipFmj^, Eq. (I4.17p . can be expressed as follows: 

7_r / (— t)"'"*"'’ r(l + air + a; 2 n) . , oi-/3i a2-/32\ /a 1 c^ 

= ^ -____^sm(7rr^+7rn^), (4.18) 

r,Ti=0 

which for f = 1, x > 0, and r = 0 (or n = 0), is identical with the series expression 
for two-sided Levy stable distribution given in, e.g. Eq. (5.8.8a) on p. 142, in [13] . 

(B) Eor 1 < aj < 2, i = 1, 2, Eq. (8.2.2.4) on p. 520 of |8], applied to Eq. (I4.16P gives 

.. cxD Ip 1 / ... j 


nr IWf / ^ ^ 'll , iiZ J-7HL—I'lp 

t '' 






M ' M m . f 

A{lp,l+j) ^ 


^mi Ip \ 
miu+lp '"'1 I 

t^yipyp i ’ 


(4.19) 


®Let us observe that, for u = M, and v = m, and thus, for a = I, and b — p, Eq. (14.1711 gives Eq. (15.111) 
of Section 5, this is the Laplace convolution of two one-sided Levy stable distributions. 
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which can be rewritten as 


H+{x,t) 


^ r!n! 

r,n=0 ^ 


r(t!r + S-')sin[™lgi-^r(„-u5)], 

(4.20) 


Example 4-1- The bi-Gaussian case. The elementary case ai = 02 = 2, Pi = (32 = 0, is 
straightforward and we include it here only for verification’s sake. Substituting Eq. 
into Eq. (14.ip and employing Eq. (3.323.2) on p. 337 of [TO], we get 

77(2, 0, 2, 0; t, x) = H e-^I+f dy = = V 2 ,oix, 2t), (4.21) 

47rr 2v27rf 

which is in agreement with Eq. (I4.20p and is presented in Fig. [1] see the curve I (red). 


Example The Gaussian-Levy case. Here ai = 2, a 2 = 1/2, /3i = 0, and P 2 
so that Eq. p4.ip reads 


77(2,0, \,-Tt,x) = 


1 ^2 7 

"I 


3 _i_ _y^ ,xy 


Vi 

dvr 

r=0 


y 2 e e 4 t 2 t dy 

(-tV4) 


r! 


f 


_3_r —ylim , 

y 2 e 4t ^ 2t dy 


1 (- 1 ) 2"+1 V+V 


= —e St 


Di 


^r!r(|+r) 21 +f’' 

^ (- 1 ) 2 -+! 


'VTt> 


2 v^ 


e St 


E 

r =0 

cx> 


(2r+ 1)! 


D 2r+1 I 
2 




1 

e St ^ " 


2 \/^ 


n=0 


n\ 


n /f3\ 4 

2 «?(-;%)> 


- 1 / 2 , 


(4.22) 

(4.23) 

(4.24) 

(4.25) 


where Di,{z) is the parabolic cylinder function [TO]. In Eq. (I4.22D we applied Eq. (9.241.2) 
on p. 1028 of [TO], and in Eq. p4.24p we changed the summation index as follows: n = 2 r+l. 
Eq. P4.25I) can be obtained from Eq. p4.19l) after using Eqs. (7.11.3.3) and (7.11.3.4) on 
p. 491 of [ 8 ]. The plot of 77(2, 0, 7, x) for 7 = 1 is illustrated in Fig. [U see the curve 

II (blue). 


Example 4-3. In Eq. (14.ip we take (7, x), j = 1,2, given in Eqs. p4.2p and (|4.5p . 
Thus, we get 


77(2,0, 



2 

exp(-|^) Re 

2 \/^ 27r 



^-iu){x-y) exp[-ifszn{i^)] 


JjJ 


77e 7°° 
2vr y_oc 
Re 7°° 
2vr y_oc 


^-iutx exp[-n^sgn(tj)] 


f 

J —c 




2\/7r7 


^-iujx exp[-::|:sgn(tj)]-Vtj’2 


do;, 


du 


dy 

(4.26) 

(4.27) 

(4.28) 
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which is in agreement with Eq. ()4.6p for a\ = 2, a 2 = 3/2, f5i = 0, and (32 = ~l/2- Thus 
we can repeat all the steps from Eq. ()4.6I) to Eq. ()4.19p which gives 


i/(2,0,|,-i;t,x) = -T5z^r(i±i + |r)cos(fj-fr) 


r=0 j=0 


r’-r- 


/l.A(3,i±i + |r)_ 

® A(6,l + j) ’ 1728t3 j ■ 


(4.29) 


Eor t = 1, the function is plotted as the curve III (green) in Eig. [TJ 



X 


Eigure 1: Multiscale densities H{ai, /3i,a2,132',t,x), for t = 1, and given values of ai, /3i, 
02, and 132 ■ In plot I (red), ai = 0:2 = 2, and f3i = [32 = 0, see Eq. (I4.2ip : in plot II 
(blue), ai = 2, /3i = 0, a 2 = ^) and /32 = — see Eq. (|4.25l) : in plot III (green), a\ = 2, 
/3i = 0, 0:2 = |, and /32 = — |, see Eq. ()4.29p . 


5 Explicit representation of the kernels of the two-scale, 
one-sided Levy generators, 0 < ct < 1 

Although this case has no direct applicability to our work on the asymptotics of the 
multifractal conservation laws in Section 3, in anticipation of our future work, in this 
section we provide explicit representations of the kernels of the two-scale, one-sided Levy 
generators with 0 < ai < 02 < 1. Here the tool is, of course, the Laplace transform. 

To further simplify our exposition we will only consider two one-parameter families, 
Vai{t,x) = ^)! and Va 2 {t,x) = ^( 02 ,- 02 , 1 )(^ 5 ^)) one-sided stable Levy den¬ 

sities whose the Laplace convolution has the form 

rx 

Vaiit,y)va2it,x-y)dy = Vaiit,x - y)va2{t,y)dy. (5.1) 

Jo 
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Functions Voij{t,x), j = 1,2 are given by the “stretched exponential” Laplace transform 
exp(—tp"-’) = C[vaj{t,x)-,p], see [25l [26l [271 [28] . Thus, the characteristic function of Eq. 
(15.11) is of the form 

exp(-tp"i = C[hai^a2{x,t);p]. (5.2) 

Eq. (j5.2p implies that, for 0 < ai = 02 < 1, 


hai,ai{t,x) = Vai{2t,x), 


(5.3) 


whereas hai,a 2 it^x) for arbitrary aj can be found by using the series representation of a 
one sided-Levy stable distribution, 


^ 00 

^{t,x) = --Yl 

TT 


{-tYT{l + ajr) . 


r=0 


ri 


^l+a^r 


3 J ' f \ 
sm(7rrajj. 


(5.4) 


Eq. ()5.4p follows from Eq. (I4.15P applied for the one-sided Levy stable distribution, and 
the formula (4) of 


Remark 2. Absolute convergence o/uQ^(f,x). Eor x > 0, the series in Eq. (15.4p 
converges absolutely and its radius of convergence is infinite. Indeed, the absolute value 
of Uc^ (x,t), for X > 0, can be estimated as follows: 


I u M ^ r(l -|- (XjU) f t 

K(Lx,)i<- 2 ; 


n=0 


ni 


X^j 


< 00 . 


Next, from the Cauchy ratio test of convergence, [29], and the Stirling formula, see Eq. 
(8.327.1) on p. 895 of [TO], the absolute convergence of the series (|5.4I) follows. The 
convergence of the series (|5.4p at x = 0 can be verified by employing the asymptotic form 
of Vaj (x, t) at X = 0: 

1 ^ 1 

^ X expf — X j, (5.5) 


where K and A are positive constants. Eq. (15. 5 p was obtained with the help of Eq. (4) 
in m- The absolute value of the right-hand side of of Eq. (15.5p obviously converges to 0, 
as X —0-|- (for any fixed t > 0, and 0 < Oj < 1). 


Eor rational aj = l/k, where k and I are positive integers, Eq. (15.41) can be expressed 
via a finite sum of the generalized hypergeometric function, see Eqs. (3), and (4), in |28j . 
for t = 1. Moreover, it turns out that, for A: < 3, it can be written down in terms of 
standard special functions, e.g., for a = 1/2, 


for a = 1/3 [3T] 


vi{t,x) = vi _i{t,x) for X > 0, 
2 2 ’ 2 


^ 3/2 / 2 i3/2\ 


(5.6) 


(5.7) 
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and for a = 2/3 [32] 


V 2 (x,t) 

3 


V TT X 2’6 



(5.8) 


where Ky{u) is a modified Bessel function of the second kind [T 0 |, and is the 

Whittaker W function [TO] , 

The substitution of Eq. (|5.4I) into Eq. (15.11) allows us to write, with help from Eq. 
(II.lip , that 


^OLl ,0.1 


/» rp OO 

it,x)= y : 

Jo 


(_t)r-+jr(l + aij)r(l + a 2 r) dy 

, rij! s.n(„au)s.n(„o,r)^ 


E 

r,j=0 


(-ty+j 


1 


r!j! r(-aij)r(-a27’) Jo 


JO 




1 {—ty~^^ r(l + aij + a2r) 

- 2^ - ^i+au-+a.r sm(7rau + vrasr). 


TT ' rli! 

r,j=0 


(5.9) 


Observe that the expressions in Eq. (15.9p are invariant with respect to the change of order 
of the parameters ai, and 02 , so that /iQi,a 2 ( 3 ^T) = ha 2 ,ai{x,t). Moreover, the first term 
in the sum in Eq. (15.9p . corresponding to j = r = 0, vanishes. The terms with indices 
j = 0 , r = 1 , and j = 1 , r = 0 , provide the ’heavy-tailed’ asymptotic behavior of the 
one-sided Levy stable distributions. Consequently, we can get see immediately that the 
asymptotic behavior of hai,a 2 it:x) is proportional to 

The first series in Eq. (j5.9jl can be expressed as a finite sum of the generalized hyper¬ 
geometric functions given in Eq. (11.91) . Indeed, let us consider, without loss of generality, 
the case of rational ai > 02 , with ai = and 02 = |, where I, k, p, and q are integers. 
In this case, Eq. (15.9p takes the form, 

hiAt,x) = - Y ^Tu--sm[ 7 r(^j + £r)] 


TT —' rli! 

r,j=0 


l+ij+^r 

X kJ q 


(_^)r+i+(n+fcpi) r[l + {{n + kpj) + ^r] 

> > > ----— - — sm\7r( in +-r) + Trlpj] 

E—/ Z—/ Z—/ r!(n -|- kpj)l -[ L q / J 


r=0 n=0 j=0 




X 




EE 


E 


"ir! r(-i„-2r)^ks 

Y[kp-i /i+i+£\ 

1 ls=0 V kp )j 


kpj 


X 


Ipj 


(5.10) 


To obtain Eq. (I5.inp we applied the Gauss-Legendre multiplication formula, and the 
Euler’s reflection formula. Now, using Eq. (II.9p we can represent the sum over j in terms 
of the generalized hypergeometric functions. Namely, for f the function hi/k,p/q{t,x) 
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can be written as follows: 


oo mi —1 

r=0 n=0 
X 


(— t)"'"’”'’ r(l + Mn + mr) 

■ sinfvrMn + vrmr) 


I r\ 


n\ r 


l+lp^mi 


^1+Mn+mr 

{Ip, 1 + 1 \ 
A(mi, 1 + n) 


1, A{lp, l + Mn + mr)_ 


(5.11) 


with K = {lpyP/m^^{t^^^/xY^, with and Mi, defined in Eq. (j4.14p . Moreover, 

it tnrns out that, after applying Eq. ( 8 .2.2.3) of [ 8 ] to the finite sum in Eq. (15.lip , we can 
write 


hi p{t,x) 

k ’ q 


X (-t)'’ /mi\M>- 

mi-lp 2^ r\ W / ^lp,mi 

{ 27 r)^- r=o ^ ^ ' 


( A(/p, 0 ) \ 

r A(mi,fr);’ 


(5.12) 


where is the Meijer G function defined in Eq. (II.7p . 

The series representing hiii^piq{t, x) given in Eq. p5.11l) converges absolutely for x > 0. 
The proof of this fact can be split into two cases: x > 0, and for x = 0. Eor x > 0, and 
fixed values of I, k, p, q, and n = 0 ,..., mi — 1 , the convergence of the series |/iz/fc,p/g(t, x)\ 
follows from the convergence of the i^ipFmi functions which, for a given n, converges to a 
constant Gn ISIE]. Thus, we get 


OO mi —1 


\hlE{t,x)\<Y, E 


^n+r 


r=0 n=0 

mi —1 


n!r! 


X 


l+Mn+mr 


sin( 7 rMn + TTmr)G„ 


^ \Gn\t"' _ 

n=0 ’ r=0 


F r(l + Mn + mr) 

^l+mr 


(5.13) 


An application of the Cauchy ratio test of convergence [29] to Eq. (I5.13P completes the 
proof of convergence for |fiz/fc,p/g(t, x)|, for x > 0 . 

Let us now show that the series representing x) converges absolutely at x = 0. 

Eor this purpose we will find the asymptotic behavior of the i+ipF^-^ function in Eq. (|5.1ip 
for «:;:$> 1, where the relation between k and x is shown below Eq. (15.111) . It follows from 
[55] . that 


l+lpPmi 


1 /" 7^\JVfn-l-TnV 1/2 —n+Mn+mr \ ™ 

m? ’ 


14) 


where, at the point x = 0, there exists an essential singularity. That gives 

1 1 r / t \^/^i - ~~r~ ^ 

hi p{t,x) I'M e ^ (—1)^” sin( 7 rMn) 

k ’ q ‘ ^ 


n M—1 ^ 

i/^mi 


n=0 


(5.15) 
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and, consequently, 


k ’ q 


mi —1 

E' 

n=0 


' 1 M-1 \ 

<mi ~^mi—lp)^ 


< i ^1+ 2(mi-ip) g- [”*l-«P+ (;4)+ 

= (5.16) 


Taking into account the fact that, for t > 0, and fixed I, k, p, and q, such that 0 < 
lfk,pjq < 1, and l/k / p/q, we can estimate Eq. (j5.16jl as follows: 

\hi p{t,x)\ < , (5.17) 

k ’ q 

where A and B are positive constants. Now , from Eq. ()5.17p it is easy to see that 
\hi/k,p/q{t,x)\ 0, for x^O. 

We will conclnde this section with several concrete examples of explicit expressions for 
bi-scale totally asymmetric Levy densities, where we can show that some of our expres¬ 
sions can be reduced to expressions in terms of more classical special functions. 


Example 5.1. The case ai = a 2 = 1/2. Substituting the corresponding Va{t,x) given in 
Eq. (j5.6p into Eq. (15.111 . we have 


hi 1 (t, x) 
2 ’ 2 


r 

4vr Jo 


dy 


P 4 y(x-y) _ 

0 [y{x - y)]^/2 


f_ r 
2vr Jo 


^/2 . 
g 4 y(x-y) 


dy 


[y{x - y)]3/2 


fOO 


dz 


TTX^ 


\Jz — \ 


t _£. 

r- 3 / 2 g =gi{x,2t). 


(5.18) 

(5.19) 


In Eq. (I5.18|l we changed the variable of integration as follows: z = 1/[1 — (1 — 2y/x)'^]. 
In Eq. (I5.I9P we also employed formula (3.362.1) on p. 344 of |10] . The same result is 
obtained from Eq. (j5.1ip where, for ai = a 2 = 1/2, we employed the formulas (7.11.3.3), 
(7.11.3.4) of [8] and (5.6.1.1) of [M]. 


Example 5.2. The case ai = 1/2 and a 2 = 1/3. Using Eqs. (15.6p and (15.7p . the Laplace 
convolution of vi(t,x), and V 2 {t,x), takes the form, 

3 3 


hi 1 (t, x) 

2 ’3 






iy . 


dy 


.El 

Ax 


[y{x - y)?/"^ 

/ + u)e-fe“ir. )<i 


(5.20) 


where u = {x — y)/y. Employing the integer representation of Ky{z) given in Eq. (8.432) 
on p. 917 in [TO] we get 
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hi 1 (t, x) 

2 ’ 3 


*2 foo g 


e 4a: / U 3 (1 + lt)6e 4a;’ 

Jo 

r -aa ^ 

,-e 4a / 3(l+t()6e 4a’ 

12v^7r2 ^ 


12v^7r2 
ai 

+o^—- 


f 


1 + u 


hJ 3 g ^ 27cca; u duJ 


du (5.21) 


w 3e “ 


l+uir 


(5.22) 


r=0 


q _ 

6 


t° \r POO 


12^3 


7r2 


r=0 


?x-3-’-(l + 7x)6+’-e-4^’^du / cj-s-’-e-'^dw 


f 


q _ 

t X 6 


-e 4a 


12v^7r2 

2“6y’^i2 

9a/ 3: 
exp(- 


E 


r=0 


_£f_\r 
273:> 

rp I 


-r(-l-r)r(-|-r)V’(-i-r;|;|i) 


-e 8^ 


X3 


CO 

%X) 


( 3^)’~^i+2r(v4;) 

r!r(|+r)r(|+r) 

x-2+2i: 


E 


■ y ~] 7 .,^^1 ( 2 x ) 3 ^ £> 1^2 

Jb 

exp(-|:) (-t) 


l+|(2+3r)(^) 

n=0 


^ (2 + 3r)! 

*2 x OO 


(5.23) 

(5.24) 

(5.25) 

(5.26) 

(5.27) 


To obtain Eq. (j5.2ip we used the series expansion of the exponential function, and there¬ 
after we changed the order of integrals and summation. We also applied the definition of 
Tricomi’s function (the confluent hyper geometric function of the second kind) tpia] 6; z) 
given in Eq. (9.211.4) on p. 1023 of [10]. In Eq. ()5.24p we employed Eq. (jl.lip and 
formula (7.11.4.8) on p. 491 of |8|, whereas in Eq. ()5.25p we used Eq. (I1.12p . In Eq. 
(I5.26P we changed the summation index from 2 -|- 3r to n. 

Note that Eq. (15.211) can be also obtained from Eq. (15.lip . Indeed, using the formulas 
(7.11.3.3) and (7.11.3.4) of |8] we get 


hi 1 (t, x) 

2’3 


E 


{-ty x"3- 


iTi 


1 + i. 


r 

Ax 


ts ’’ r{- 5 ) 

„-t^/(8a:) °° ( f\r 


E 

r=0 


i-t) 


. r 3 

’■+^ x"3-2 


’’ r(-j-a 


y/^: 


TTX “ rl 
r=0 


iFi 



(5.28) 
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Example 5.3. The case ai = 2/3, a 2 = 1/3. Here, the Laplace convolution defined in Eq. 
m is equal to 


hi 2 (t, x) 

3’3 


3 

t2 




\/37r2 Jo 

2\/3 




/■°° ex 

Ji 


exp(- 


2t^ 

ThP- 


U 


(5.29) 


,2*3/2 


1/4 


(V^- 1)3 


9 _£ 

3 

127ri 


/■°° ex 


i / ^ ,5 

(-/U — 1)3 


i( 

^ 9 ’ R 


4/3 


t^x 3 /■°° exp(-^u) 


127ri 


2^-1 °o 


e{^/u- 1)3 2 6 


t X 3 


r°° ex 

Ji 

°° ( \r fc 

^—n * *^0 



4 


y/u 


UJ 3e 27 x<.; 


E 

r=0 


27xui yiI-1- 


^ dcj 


r! 


4 

603 


du 


a; 3 


doj 


1277 2 

^ (-ife)' 


E 


1277 2 


t\ 


r(— 1 — r)Int(x,t). 


(5.30) 


Above, the integral representation of K^{z), see Eq. (8.432) on p. 917 in [TO], was used. 
The auxiliary function Int(x,t) is defined as follows: 


Int(t,x) = / u 6 ’'' 2 (y^—1) 3 + 2 e 27x2“ lEi 1 du 



(5.31) 


where Eq. (2.19.1.4) on p. 201 of [8] was employed. Substituting Eq. (15.311) into Eq. 
(I5.30D . using the Euler’s reflection formula, the Gauss-Legendre multiplication formula, 
and changing the summation index as follows n = 2 + 3r, we get 


hi lit,x) = ^ 


71=0 


n\ r(-f)'^' 


A _l_ n 1 I 71 
2' fi ’ “I” fi 


6 • 

1 2 
3’ 3 



oo 


+ E 

71=0 


(-t)l+’3 

n! 



n 

3 



2^2 


5 I n 

6 6 ’ 


4 
3 
2 4 
3’ 3 


+ - 
^ 6 . 



oo 


+ E 

71=0 


o I —l — H 

(-t)7+" X 3 3 
2n! r(-J-S)"^ 


f 1 + n 5 + 

( 6 ^ 6 ’ 3 ^ 
I 4 5 

\ 3’ 3 


n 

6 . 



(5.32) 
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Finally, we would like to point out that 2 -^ 2 ’s cannot be expressed in terms of any 
other special functions. And, of course, the explicit form of hi 2 (x,t) can be obtained 

3’3 

from Eq. (15.111) and Eq. (15.321) . In Eig. 5.2 we show plots of the densities discussed in 
Examples 5.1-3. They were efficiently obtained via the symbolic manipulation platform 
Mathematica 



X 


Eigure 2: The plots of hai,a 2 {x,i)- Plot I (red): cci = 0:2 = see Eq. (I5.18p i: plot II 
(blue); «! = 02 = |, see Eq. (15.281) : and plot III (green): cci = |, 0:2 = see Eq. 

( 1021 ) . 


6 Conclusions and comments; the Stieltjes moment problem 

The goal of this work was to show how the asymptotic behavior of solutions of certain 
nonlinear conservation laws can be explicitly described in terms of some special functions 
which make the efficient computation of the related probability densities possible. In our 
future work we plan to expand this work to other types of asymptotic problems for other 
nonlinear evolution equations discussed in, e.g., m, and other papers referenced therein. 

Here, we want to conclude with an example of how the results obtained in Section 5 
can be applied to the solution of the classical Stieltjes moment problem for special cases 
considered in Examples 5.1-3. Recall, that the Stieltjes moment problem [35l [36] can be 
formulated as follows: Find a positive function W(x) which satishes the infinite set of 
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equations, 

x^W{x)dx = p{n), n = 0,1,2,..., (6.1) 

for a given moment sequence {p{n)}. 

From the practical point of view, it plays a special role in probability theory where the 
moment estimators can be usually conveniently obtained, and the issue is whether they 
determine the relevant probability distribution. In general, the solution to the Stieltjes 
moment problem can be unique or non-unique. The examples of non-unique solutions can 
be found, e.g., in [36j . Below, we quote one uniqueness criterion, and one non-uniqueness 
criterion to be used below. 

Cl Carleman uniqueness criterion, (see, .e.g., [35]). If S = Yl'^=i[p{n)] 2 n = oo, then 
W{x) is uniquely determined by Eq. (16.11) . 

C2 Carleman non-uniqueness criterion, see, .e.g., |37l[38]. If Yl^=i[p {'^)]~< oo, and 
if there exists x' > 0 such that, for all x > x', W{x) > 0, and '4^{y) = — ln[iy(e^)] is 
convex in y' > 0, where y' = In(x'), then W{x) is non-unique. 



Let us begin with finding the /i-th moments of hai,a 2 ixjt)' 


Pai,a2ip) = / xf^hai,a2it,x)dx, pG 
Jo 


( 6 . 2 ) 


Using Eq. (|5.12j) . changing the variable x to y = x applying Eq. (2.24.2.1) in [8], and 
the Gauss-Legendre multiplication formula (see Eq. (|1.12p l. we arrive at 


Pa\,a2{p) ~ 


1 r 

Ip Jo 


— 1 
y ‘p h 


ai ,a2 


{y 'P,t)dy 


M ^ 

r=0 


t>0. 


r\ 


(6.3) 

(6.4) 


The absolute convergence of the series in Eq. (16.411 is easy to check through the Cauchy 
ratio test. Eq. dsai) implies that the /ith moment is finite for a given r, and p < mr, 
whereas it is infinite for p > mr. The moment Pai,a 2 {p) is equal to one for p = 0. The p- 
th moment of /iQ,^^Q, 2 (x,t), for p = —Ipn, n = 0,1, 2,..., gives rise to the Stieltjes moment 
problem for which, from the comparison of Eq. (16.31) with Eq. (j6.1h . W{x) and p{n) are 
of the form, 

1 - 1 -— — — 

Wa^^a 2 {t-,x) = —X ''Phai,a 2 {'t,x ‘p), and p(n) = (-/pn). (6.5) 

Now, we can employ the above criteria Cl and C2 to take a closer look at the Stieltjes 

moment problem related to the Examples 5.1-3. 

In Examples 5.1-2 functions Wa^^a 2 {x) are unique. Indeed, for ai = a 2 = 1/2, we 

have p{n) = 2(2t)“^’^r(2n)/r(n), which gives S = 2tY^^=i[‘J^{‘Jn)/T{n)]~^. Using the 

Stirling formula for the gamma function, Eq. (8.327.1) on page 895 of [l0|, we get that the 

series 5 Ri t-y/e/(2\/2) n the latter series is divergent. Thus, Cl leads to the 

unique function VFi i(t,x) = t/(-GTO) exp(—xt^). Similar considerations are valid of the 
2 ’ 2 
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case of «! = 1/2 and a2 = 1/3 which leads to the uniqueness of the appropriate function 
W{t, x) in that case as well. 

However, Example 5.3 leads to a non-unique Stieltjes moment problem. Indeed, 
for ai = 2/3, and a 2 = 1/3, we have p{n) = 3t“®”'r(6n)/r(2n), and the series S = 
3“^[r(6n)/r(2n)]“^ « t^e^/108is convergent (we used the Stir¬ 
ling formula Eq. (8.327.1) on page 895 of [10] here). The second condition in C2 requires 
verification of the sign of the second derivative of ihi 2 (t,x) = — InlWi 2(t,e^)l, 

3’3 3’3 


(P' 4’1 2 (t, x) 

3’3 

dx"^ 


d'^ ln[ITi 2(t, e^)] 




d2 


= —uhl lit,u) —hi 2 (t,u)hi 2 (t,u) 

3 ’ 3 du 3^3 3 ’ 3 

1 2 


uhl 2{t,u)—^ + 

3 ’ 3 du^ 


d , / N 

— hi 2 (t, u) 

du 3’3 


( 6 . 6 ) 


for u = exp(—x/2). If Eq. (|6.6I) is positive then i/jiy) is convex and lEi 2 (t,x) is non- 
unique. The positivity of Eq. (|6.6I) with hi 2 {t, u) given in Eq. (15.321) can be established 
analytically but for the purpose of this commentary section we are just showing the plot 
of Eq. (16.61) in Fig. [3l 



;c 


Figure 3; The plots of i{t,x) for different values of t. Plot I (red): t = 0.8, plot II 

^ 3^3 

(blue): t = 1, and plot III (green): t = 1.2. 
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